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^ ; ABSTRACT 

H , I discuss the possibility that accreting, supermassive black hole (SMBH) binaries with 

' sub-parsec separations produce luminous, periodically recurring outbursts that inter- 

rupt periods of relative quiescence. This hypothesis is motivated by two characteristics 
found in simulations of binaries embedded in prograde accretion discs: (i) the formation 
of a central, low-density cavity, and (ii) the leakage of gas into this cavity, occurring 
once per orbit via discrete streams on nearly radial trajectories. The first feature will 
diminish the emergent optical/UV flux of the system relative to active galactic nuclei 
I (AGN) powered by single SMBHs, while the second can trigger periodic fluctuations 

r ) • in the flux. I propose a simple toy model in which a leaked stream crosses its own orbit 

' and shocks, converting its bulk kinetic energy to heat. The result is a hot, optically 

[ thick flow that is quickly accreted and produces a flare with an AGN-like spectrum 

Oh, that peaks in the UV and ranges from the optical to the soft X-ray. Due to the pre- 

Q ' ceding quiescence, such a flare could plausibly be mistaken for the tidal disruption of 

i—{ , a star. For typical binary periods of years to decades, the event rate in an individual 

tj^ • system can be much higher than that predicted for stellar tidal disruptions but infre- 

d I quent enough to hinder tests of periodicity. The flares proposed here can be produced 

by very massive (> 10^ M©) SMBHs that are thought not to tidally disrupt solar-type 
stars. They could be discovered serendipitously in the future by observatories such as 
^ I LSST or eROSITA. I apply the model to the active galaxy OJ 287, whose production 

Q>^ . of periodic optical flares has long fueled speculation that it hosts a SMBH binary. 
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INTRODUCTION If identified, observing such systems can unveil invalu- 

able clues about the physics of accretion processes in ac- 



^ Given that massive galaxies assemble through hierarchical , i ■ / a^-iatn j j_i i j_- r i • 

' , , . n , , tive galactic nuclei (AGN) and the co-evolution oi galaxies 

i. _j 1 mergers and that virtually all ot them host a supermas- , ,i ■ i nuT-nTT a ■ h , . i- ■ -i -i 

jic' . ,,,,,,„, ,r,TT\ . , . , • / 1 ■ and their nuclear bMBHs. An especially tantalizing possibil- 

' sive black hole (bJVlBH) m their nuclei (see, e.g., the review r i ■ i ■ ■ ■■ .i 

?H ■ , „ V, 1 n r , \ ity IS that OI observing near-merger binaries concomitantly 

rrt , by Ferrarese & Ford 2005, and references therein), it is m- 1^1,^11^ ^- j i- 1 

^ , , „ „ , , through both electromagnetic and gravitational- wave signa- 

evitable that they spend a fraction 01 their time hosting two , at , 1 11 1 i,- . i- -1 

tures. Not only would such multi-messenger studies provide 
or more SMBHs. Such a conhguration will eventually re- , . , ■ 1 , . 1 . ... r , 1 1 

, . , r ■ r , 7o»«"i-,TT , • , ,. unprecedented independent determination ot the masses and 

suit m the formation ot a bound SMBH binary, whose orbit . . , , . , .^-iat ■ ^^ 11 1 ^ 

spins 01 the central AGN engines, they can also be used to 
becomes increasingly compact as it loses its orbital angular , • ■ , ■ ^ /tt 1 tt 1 nr^r^r 

, „ , , , . , ^ . . probe the cosmic expansion history (Holz & Hughes 2005; 

momentum and energy-(I) first through dynaniical friction ^^^^.^ 2006; Bloom et al. 2009; see Schnittman 2011 

with the stellar and gaseous background, (II) then through „ ■ r 1 1 ^ ^- ^ ^ \ 

" , ' , , tor an overview or proposed electromagnetic counterparts), 

three-body interactions with stars and/or tidal interactions 

with a gaseous accretion disc, and (III) finally, at separa- SMBH binaries, however, have thus far mostly eluded 

tions well below a parsec, being driven to merger through discovery. Although there are now numerous examples of 

the emission of gravitational waves (Begelman, Blandford & dual and even triple AGN at kpc-scale separations (e.g., 

Rees 1980; Yu 2002; Haiman, Kocsis & Menou 2009; Colpi Komossa et al. 2003; Bianchi et al. 2008; Comerford et al. 

& Dotti 2011). 2009; Green et al. 2010; Liu et al. 2010; Barth et al. 2008; 

Liu, Shen & Strauss 2011), there is only one unambigu- 
ous example of a proper — that is, gravitationally bound — 

* E-mail: taka@mpa-garching.mpg.de binary, which has a projected separation of 6 pc (Rodriguez 
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et al. 2006). In searches for periodic variability in AGN, a bi- 
nary explanation has been ruled out for dozens of candidates 
(Halpern & Filippenko 1988; Eracleous et al. 1997; Halpern 
& Eracleous 2000; Eracleous & Halpern 2003; Chornock 
et al. 2009). While there are several examples of AGN with 
double-peaked broad lines (Gaskell 1983; Peterson, Korista 
& Cota 1987; Gaskell 1996; Boroson & Lauer 2009), shifted 
lines (Bogdanovic, Eracleous & Sigurdsson 2009; Dotti et al. 
2009), or suggestive radio lobe morphology (Roos, Kaastra 
& Hummel 1993; Gopal-Krishna, Biermann & Wiita 2003; 
Liu, Wu & Cao 2003), such features have alternate explana- 
tions that do not require the presence of a binary (see the 
review by Komossa 2006). The blazar OJ 287 (a.k.a. EGO 
0851-1-202), which we will revisit in §3.3, has been suggested 
to host a binary SMBH because it produces pairs of optical 
flares at a regular interval of approximately 12 years (Sillan- 
paa et al. 1988; Valtonen et al. 2008). 

The dearth of observational evidence is not surprising, 
for several reasons. First, the detection of SMBH binaries 
must lean heavily on indirect observations of kinematic sig- 
natures over long periods of time, since resolving images of 
sub-parsec binaries is intractable in all wavelengths but ra- 
dio, where it is feasible only for relatively nearby sources. 
Second, because gravitational wave emission and accretion 
disc torquing accelerate the orbital evolution of compact bi- 
naries, the most compact binaries — whose kinematic signa- 
tures would be the easiest to detect — should also be the 
rarest (e.g. Haiman, Kocsis & Menou 2009, and references 
therein). Third, the emission features of an accreting SMBH 
binary may deviate significantly from those of AGN host- 
ing a single SMBH, potentially causing such systems to be 
missed in AGN surveys. For example, many theoretical stud- 
ies suggest that if the binary's SMBH masses are similar and 
its orbit is prograde with its accretion disc, then the binary's 
tidal torques can suppress the accretion rate in its vicinity 
and form a low-density, central, circumbinary cavity in the 
disc (e.g., Artymowicz & Lubow 1994; Hayasaki, Mineshige 
& Sudou 2007; MacFadyon & Milosavljevic 2008; Farris et al. 
2012). Such a cavity would result in reduced emission in soft 
X-ray, UV and even optical frequencies compared to typical 
AGN (Tanaka, Menou & Haiman 2012; Gultekin & Miller 
2012). 

Numerical sirrmlations of the cavity agree that the sup- 
pression of the central accretion rate is not absolute. Once 
per orbital period, gas leaks into the cavity in discrete, elon- 
gated streams (e.g., Artymowicz & Lubow 1996; Cuadra 
et al. 2009; D'Orazio, Haiman & MacFadyen 2012), a por- 
tion of which will accrete onto one or both SMBHs. If the 
accretion timescale of an individual stream is long, then the 
leaked gas will form a small accretion disc around one or 
both SMBHs (Hayasaki, Mineshige & Ho 2008). 

In this work, I consider the possibility that the accre- 
tion timescale of the gas leaked into the cavity is short, due 
to the conversion of the stream's excess orbital energy to 
heat. In this scenario, the heated gas fuels a hot, optically 
thick accretion flow that results in a luminous flare. The 
idea that leakage of gas into the cavity might cause variable 
AGN emission is not new; the strongly modulated accretion 
rates onto either SMBH found in simulations have long been 
invoked to explain objects such as OJ287 (Artymowicz & 
Lubow 1996; Hayasaki, Saito & Mineshige 2012). The main 
goal of this study is to point out that such accretion events 



are likely to exhibit transient optical/UV flares that (i) can 
occur in galactic rmclei that appear quiescent and (ii) should 
be discovered, potentially in large numbers, by wide-angle, 
high-cadence surveys such as LSST^ or eROSITA^ . If iden- 
tified, very compact accreting SMBH binaries will present 
exciting targets of study in both gravitational and electro- 
magnetic waves (cf. Haiman et al. 2009; Tanaka, Menou & 
Haiman 2012; Sesana et al. 2012). 

I present a toy model to demonstrate that for plau- 
sible estimates of the stream mass and viscosity, the gas 
accretes to the nearby SMBH and becomes depleted and 
optically thin within an orbital period, i.e. before the next 
stream is deposited. In this case, such flares would be tran- 
sient events that recur once or twice every orbital period, 
depending on whether one or two SMBHs produce a flare. I 
show that if the sound speed in the accreting gas is propor- 
tional to the Keplerian orbital velocity, then the accretion 
rate onto the SMBH will decay as roughly M oc t~*^^ until 
the flow becomes optically thin. Such flares could plausibly 
resemble stellar tidal disruption events (TDEs) whose ac- 
cretion rates have a canonical power-law decay A4 oc t~^^^ 
(Rees 1988; Phinney 1989). Unlike tidal disruptions, how- 
ever, these flares can occur in pairs (if both SMBHs ac- 
crete), and in more meissive (> 10* Mq) SMBHs that would 
not disrupt Sun-like stars. 

This paper is organized as follows. In §2, I review the 
overall geometry of the accretion flow onto the SMBH bi- 
nary, namely (i) the central cavity in a prograde circumbi- 
nary disc and the corresponding reduction in the emergent 
spectrum, and (ii) the leakage of discrete streams into the 
cavity. A semianalytic toy model for the accretion of the 
stream and the subsequent time-dependent, spectral emis- 
sion are presented in §3. I apply this model to produce mock 
light-curves and spectrum snapshots of flares for specific ex- 
amples: a 10'^ Mq SMBH binary (§3.2), and a 10® M© bi- 
nary whose orbital parameters are motivated by the blazar 
OJ 287 (§3.3). Upper limits to the frequency and popula- 
tion of flaring SMBH binaries are estimated in §4, demon- 
strating a proof-of-concept that such flares might be quite 
common. Observational implications, including the similar- 
ity to TDEs, as well as theoretical caveats are given in §5. I 
conclude by summarizing the findings in §6. 



2 ACCRETION FLOW GEOMETRY 

Let us consider a prograde accretion disc around a SMBH 
binary of mass M,,, semimajor axis a and orbital period 
P. The Lindblad resonance clears the gas around the bi- 
nary's orbit, and inhibits the infall of gas inward of a radius 
2a, excavating an annular gap (Artymowicz & Lubow 
1994; but such a gap may close — see Kocsis, Haiman & 
Loeb 2012b,a). For binaries with comparable BH masses, 
the gas can cross inward of « 2a only once per orbit on 
nearly radial trajectories, at a somewhat suppressed time- 
averaged rate — typically ~ 10 — 100% of the accretion rate 



^ Large Synoptic Survey Telescope 
http : / /www. Isst . org/lsst/ 

extended ROentgen Survey with an Imaging Telescope Array 
http : //http : //mmm . mpe . mpg . de/eROSITA 
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at larger radii (e.g. Ochi, Sugimoto & Hanawa 2005; Mac- 
Fadyen & Milosavljcvic 2008; Cuadra ct al. 2009; D'Orazio, 
Haiman & MacFadyen 2012). It is convenient to think of the 
accretion flow in terms of two dynamically distinct compo- 
nents separated by the gap^: an outer, geometrically thin 
and optically thick circumbinary disc at > 2a, and an in- 
ner, low-density region a.t R ^ 2a. In this section, I present 
a rough physical description of each component. 

2.1 The outer disc 

The outer disc, to a first approximation, may be thought of 

as a standard accretion disc around a single SMBH but trun- 
cated sharply inside R ~ 2a. I will ignore here the certain 
aspects of the cavity-disc boundary, such as its non-circular 
shape, precession, and perturbations in pressure and tem- 
perature (MacFadyen & Milosavljcvic 2008). 

Because it is "missing" the innermost regions that pro- 
duce high-energy thermal photons in an AGN disc, the cir- 
cumbinary disc is intrinsically dim at optical and UV fre- 
quencies when compared to a disc around a single SMBH 
with the same mass. Under the simple assumption that disc 
is a superposition of annuli that emit as isotropic graybod- 
ies, the spectral energy distribution (SED) L„ is given by 

L^ = 2 1 2nR F^{Tp) dR, (1) 

where F^Tp) = 7r-^S,(rp) (2) 

is the emergent flux per disc face, is the ratio of the ab- 
sorption opacity to the total opacity (see, e.g., Blaes 2004), 
and is the Planck function. 

The effective temperature Tp{R) of the thermalization 
photosphere is given by 

E{R,T,)asBT^{R) = §;^(^l-^l^^Mnl. (3) 

Above, -Risco is the radius of innermost stable circular or- 
bit, and the dimensionless factor H(_R, Tp) (Tanaka & Menou 
2010) quantifies the deviation of the emission from a black- 
body; at low temperatures (large disc radii), H « 1. The 
Keplerian angular velocity is Qk, and M is the mass accre- 
tion rate, which I will parameterize as 

M = rh^, (4) 

with m a dimensionless parameter quantifying M in 
terms of the Eddington luminosity I/Edd(M) = 1.25 x 
10'^** M/M0ergs^^ and a radiative efficiency e. I adopt 
e = 0.1 throughout this paper. 

In a steady-state disc around a single SMBH, the stan- 
dard assumption is to take M to be constant with radius. 
This is not the case for the circumbinary disc; outside the 
cavity, the right-hand side of equation (3) becomes modi- 
fied due to the deformation of the surface density profile 
by the binary torques (e.g. Ivanov, Papaloizou & Polnarev 

^ These components are not completely independent of each 
other. The gas that leaks into the cavity can return to and interact 
with the disc-cavity boundary; D'Oraaio, Haiman & MacFadyen 
2012. 
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Figure 1. Model accretion disc SEDs with m = 1, e = 0.1. The 
spectrum that extends to the highest frequencies in each panel 
shows the emission from a graybody disc extending inward to 
R = 3GM../c2, as in a single SMBH engine. The less energetic 

spectra correspond to discs with the same central engine mass 
and accretion rate, but truncated by a cavity inside R = 2a. The 
anomalous optical, UV and soft X-ray emission may cause some 
accreting SMBH binaries to escape AGN surveys. 

1999; Liu & Shapiro 2010). I should note that because the 

central cavity strongly suppresses the disc's luminous out- 
put, the parameter rh can greatly exceed unity without the 
bolometric luminosity exceeding the Eddington limit; that 
is, the surface density can be substantially higher than in 
AGN accretion discs around single SMBHs. 

For a disc truncated inside R~2a, the emergent spec- 
trum i/Lv peaks at a frequency t'peak ~ (2 — 3) x fcB2p(2o)//i, 
then drops off steeply at higher frequencies. Assuming H w 1 
and a S> i?isco, this peak frequency evaluates to 

lo" (life) (k^^)""«^^ 

Figure 1 shows the model graybody spectra for accretion 

discs around SMBHs with masses of M = 10*^ M© (top 
panel) and 10* M© (bottom) and m — 1. The curves show, 
in decreasing order of peak photon energy, a graybody disc 
that extends inward to Risco (taken here to be 3GM,,/c^) 
and those that are truncated at 2o, with o corresponding to 
binary periods P — O.lyr, lyr and lOyr. 

The most prominent consequence of the central cavity 
is that for a wide range of binary separations, the accre- 
tion disc is dim at UV frequencies, unlike typical AGN (see 
Tanaka, Menou & Haiman 2012; Gultekin & Miller 2012). 
The majority of SMBH binaries arc expected to have long or- 
bital periods of years to centuries (Yu 2002; Haiman, Kocsis 
& Menou 2009), so the thermal disc emission will typically 
be cut off at UV, optical, and perhaps even infrared wave- 
lengths. Accreting SMBH binaries may therefore be missed 
by AGN detection methods that rely on UV and optical 
luminosities or colors. If soft X-rays in AGN are also pro- 
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duced in the central regions of an accretion disc (e.g., Done 
et al. 2012), then some binaries could also elude X-ray AGN 
surveys. The low central gas densities may also result in 
weakened broad-line emission relative to typical AGN, and 
contribute to a hard X-ray corona. 

2.2 Central cavity and leaked streams 

The SEDs presented in Figure 1 represent only a part of the 
total emission. As stated above, the outer disc leaks gas peri- 
odically into the cavity, and to predict the thermal emission 
from the whole system we must understand what happens to 
this leaked gas. Numerical simulations show that the leak- 
age takes place in discrete streams that enter the cavity once 
per binary period on nearly radial free-fall trajectories and 
preferentially accretes onto the secondary (Artymowicz & 
Lubow 1996; Cuadra et al. 2009). 

How far from either SMBH is the stream pericenter? In 
a simulation of accretion onto a SMBH binary with a 1 : 3 
mass ratio, Cuadra et al. (2009) found the accretion across 
a radius of 0.1a from either SMBH to be highly punctuated 
and periodic, indicating that pericenter lies inside this ra- 
dius. Sesana ct al. (2012) resolve in their simulations of a 
binary with a 0.35 : 1 mass ratio small accretion discs with 
radii ^ 0.005a around both SMBHs, which suggests a similar 
upper limit to the pericenter radius. (That is, the inner discs 
would be larger if the pericenters of the streams were farther 
from the SMBHs. Note that 0.1a is smaller than the canon- 
ical Hill stability radius, Rmu ^ a(l - e)[Af./(3Af..)]^/^, 
where M, is the mass of either SMBH and Af.. is the to- 
tal binary mass.) Similarly, the simulations of Hayasaki, Mi- 
neshige & Ho (2008) for a 1 : 1 system — under the assump- 
tion that the streams cool efficiently and the entire accretion 
flow onto the binary is isothermal — form small, dense discs 
with radii ~ 0.1a. 

If the streams accrete on timescales longer than the bi- 
nary period, then over time the leaked gas will sustain small, 
optically thick accretion discs around one or both individ- 
ual SMBHs (e.g., Hayasaki, Mineshige & Ho 2008). In this 
scenario, the inner disc(s) would emit in the UV, but more 
weakly than in single-BH AGN due to the lower gas content. 
The binary could also exhibit periodic optical/UV flares as 
the incoming periodic streams interact with the pre-existing 
inner disc(s) (Tanaka, Menou & Haiman 2012). 

An alternate possibility, which I explore further in this 
manuscript, is that the leaked gas is consumed by one or 
both SMBH(s) on a timcscale faster than the binary's or- 
bital period. In this case, no persistent discs are present 
inside the cavity, and the streams instead fuel transient ac- 
cretion flows that produce UV-optical flares once per binary 
orbit. In between streams, the accretion flow would exhibit 
low intrinsic UV flux and unusual spectral shape like the 
examples shown in Figure 1, and the system would appear 
as an odd AGN or even masquerade as a quiescent nucleus. 
A key consideration is the large amount of kinetic energy 
that the streams have at pericenter; the gravitational po- 
tential difference between the cavity edge (~ GM„/{2a)) 
and pericenter ( > GAf,/(0.1a)) is considerable. 

The low-density cavity would be optically thin, and 
whatever heat is generated there would be unable to cool 
efficiently. Thus, in between the periodic leakage of the 
streams, the central cavity should resemble the underlu- 



minous, ADAF (advection dominated accretion flow) state 
(Ichimaru 1977; Narayan & Yi 1994, 1995). It is perhaps 
interesting to note that the conditions at the boundary be- 
tween the gap and the circumbinary disc — such as steep 
gradients in temperature, pressure and radial velocity, as 
well as super-Keplerian orbital velocities near the bound- 
ary (MacFadycn & Milosavljcvic 2008) — satisfy the theo- 
retical requirements for consistent solutions with an "outer 
disc, inner ADAF" geometry (Honma 1996; Liu et al. 1999; 
Manmoto &: Kato 2000). Despite producing virtually zero 
emission in the optical to soft-X-ray frequencies, the opti- 
cally thin accretion flow inside the cavity may help generate 
outflows (Narayan & Yi 1995; Blandford & Begelman 1999; 
Meier 2001; Tanaka & Menou 2006; Narayan & McCUntock 
2008), which may be observable as radio sources or as blazars 
(e.g., Falcke & Biermann 1999; Cao 2002). 



3 MODELING THE PERIODIC FLARES 
3.1 Viscous spreading 

I now turn to modeling the accretion of the stream and 
the resultant emission. The stream is tidally elongated as 
it falls, and upon passing the point of nearest approach to 
a SMBH, the leading end will collide with the rest of the 
stream, shock-heating the gas; additional heating may occur 
due to tidal compression, as well as the enhanced velocity 
gradients and viscous shear at pericenter. The expected re- 
sult is that a moderate fraction of the stream's initial bulk 
kinetic energy is converted to radiation: 

^ (l - ^1 - eias) P^K ~ aradT*, or (6) 

^pwK ~ 3pcs. (7) 

Above, the quantity vk is the circular Keplerian velocity at 

the radius where the shock is formed, and the left-hand side 
of equation (7) is the kinetic energy in excess of that for 
a circular orbit; egas ;^ 1 is the orbital eccentricity of a gas 
element in the stream before the collision, p is the gas volume 
density, T is the gas temperature and Cs is the isothermal 
sound speed. As noted above, the amount of kinetic energy 
that is available to heat the accretion flow is very large. 

What is the mass Afflare that is accreted in this fashion? 
As stated above, the time-averaged mass flux into the cavity 
is of order ~ 10—100% of the accretion rate in a disc around 
a single SMBH in simulated examples, so the mass that en- 
ters the cavity per leakage event is ~ (0.1 — l)MR>2a P- 
Presumably, only a fraction of this mass circularizes via 
shock heating and accretes onto the SMBH, with the rest 
being flung back out into the circumbinary disc (D'Orazio, 
Haiman & MacFadyen 2012, and contribute to a pile-up 
of gas just outside the cavity as discussed by Ivanov, Pa- 
paloizou & Polnarev 1999, Liu & Shapiro 2010 and others) 
or contributing to an outflow. I parameterize the mass as 

Aiflare = fMR>2a P 
P 

= frhR>2ai — M,, 

ISal 
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where / ~ 0.1 quantifies the mass that ends up in the accre- 
tion fiow, and isai ~ 45 Myr is the Salpctcr time. Whatever 
leaked mass is not accreted immediately may be re-ejected 
into the circumbinary disc (see earlier footnote), fall back 
onto either SMBH at a later time, interact with and repro- 
cess the photons produced by the flare, contribute to the 
ADAF/corona within the cavity, or ejected via thermally, 
radiatively or magnetically driven outflows. While the un- 
accreted gas will plausibly have observable features, the flare 
must be the energetically dominant event, at least in the UV 
and soft X-ray frequencies, because it is the most signiflcant 
component that can tap into the deepest parts of the bi- 
nary's potential well. 

Given the range of dynamic scales of the problem and 
the number of relevant physical processes involved — the ax- 
ial asymmetry and fully three-dimensional aspect of the ac- 
cretion event, steep thermodynamic and density gradients, 
the importance of photon forcing and photon trapping, a 
magnetohydrodynamic treatment of the viscosity mecha- 
nism, to name a few— it will be a computational challenge 
to provide a detailed theoretical model of the flare. Instead, 
I present a toy model, in which the fluid elements have Ke- 
plerian circular orbits"' and the dynamics can be modeled as 
a geometrically thin disc by integrating the relevant quan- 
tities in the vertical direction. This approach is admittedly 
crude, but serves a heuristic role by providing (i) a proof- 
of-concept that a fairly conservative viscosity prescription 
can drive a rapid, optically thick accretion event, and (ii) 
an analytic foundation within which future, more detailed 
studies can be framed. 

The shock-heated and gravitationally bound gas will 
spread as turbulent viscosity transports angular momentum 
and dissipates orbital energy. The kinematic viscosity f is 
the product of the gas sound speed and the size of the largest 
turbulent eddies. The latter cannot be larger than the disc 
scale height H = Cs/H, so I parameterize the viscosity in the 
familiar fashion, v = aHca = ac^/O,, with a < 1 (Shakura 
& Sunyaev 1973). Empirical estimates of a in AGN typi- 
cally give a ~ 0.01 — 0.3, while outbursts from accreting, 
stellar-mass compact binaries give a ~ 0.3 (e.g., Dubus, 
Hameury & Lasota 2001; Starling et al. 2004; King, Pringle 
& Livio 2007). Simulations of the magnetorotational insta- 
bility suggest similar values (Pessah, Chan & Psaltis 2007). 
Dissipations of supersonic shocks, such as the ones consid- 
ered here, or other instabilities could drive macroscopic tur- 
bulence with a ~ 1 (e.g., Carciofl et al. 2012). In magneto- 
hydrodynamic simulations of circumbinary discs, Shi et al. 
(2012) and Noble et al. (2012) flnd Reynolds and Maxwell 
stresses equivalent to a > 1 near the cavity edge. Using 
equation (7), I prescribe 
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a' 

6 CIk 
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(9) 



The above expression for v is tied to the gravitational poten- 
tial of the nearby SMBH but not explicitly linked with the 
thermal content of the gas; it is equivalent to the standard 
ShakurarSunyaev prescription only as long as Cg ~ i^k/S. I 
denote the viscosity parameter as a' for distinction. I will 

* The present analysis should not be significantly afEected if the 
bound gas forms an eccentric accretion disc; see Syer & Clarke 
(1992); Lyubarskij, Postnov & Prokhorov (1994); Ogilvie (2001). 



calculate the viscous evolution using equation (9) instead of 
the standard prescription f = ac^/flK- I discuss the utility 
and limitations of this approach in §3.2 below. 

I use the standard equation for Keplerian, viscous thin 
discs (e.g., Pringle 1981), which gives the evolution of the 
vertically integrated density E: 



27vR^^E{R,t) 



_d_ 
dR 



(10) 



The viscosity form suggested in equation (9), p (x i?'''^ has 
well-documented exact solutions for arbitrary initial proflles 
of S (Liist 1952; Lynden-BeU & Pringle 1974; Pringle 1991; 
Tanaka 2011), which I will employ throughout the rest of 
this paper. 

Initially, our accretion flow will be moderately geomet- 
rically thick, with H^/R^ ^ cl/v^ ~ 1/6. The pressure is 
dominated by radiation, and the main source of opacity is 
electron scattering. 

The infall timescale of the gas is 



U{R) = 



R_ 

Vr 



4_ 
a' 



1 + 2- 



;in(:/S) 



dlnR 



(11) 



For an initial mass distribution initially concentrated at a 
radius Ro, 



U{Ro) 



< 



3/2 



(12) 



That is, the accretion flow will be depleted on a timescale 
shorter than the binary orbital period (cf. Sesana et al. 2012) 
if it is initially deposited at a radius 

a 



Ro — 



M, 



1.0 



0.4 



a 

or 



2/3 



Rb 



(13) 



Above, I set e ~ 0.6 as a reference value for the binary 
eccentricity (Roedig et al. 2011). As discussed above, this 
criterion appears to be amply satisfied by the orbital tra- 
jectories of leaked streams in numerical studies. The actual 
accretion timescale might be even shorter if thermal or vis- 
cous instabilities drive up the viscosity (see §3.2 below). 

Note that the above relations imply a simple scaling 
of the viscous evolution with the system parameters. For 
any given values of the binary mass ratio and a' , and a 
set surface density distribution with respect to the dimen- 
sionless variable R/a, the viscous evolution with respect to 
the variable t/P is identical regardless of the values of M,, 
and a. This invariance will break down in regions where the 
accretion flow becomes optically thin, where geometrically 
thin, radiatively efficient solutions cannot exist. I compute 
T = Oites^, where 6 = 0.2 is a porosity factor (Turner 2004; 
Milosavljevic & Pliinney 2005), and neglect the emission 
from the accretion flow where r < 1. Because the accre- 
tion flow spreads so that the surface density at large radii 
is lower than near the center (see Tanaka 2011), the outer 
regions of the disc become optically thin first. 

Let us now apply this model to a specific example. I 
consider a binary of total mass M,, = IO^Mq, with a period 
of P = 30 yr (a = 0.010 pc). I set mR>2a = 1 and / = 0.05, 
setting the mass in the accretion flow as A/fiarc = 0.33 Mq. I 
assume accretion onto the secondary SMBH with M,/M,, = 
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1/10 and Jiisco = SGM,/c^. In the top panel of Figure 2, I 
plot the mass flux 



0.001 



M{R,t) = 3nvT. 



1 + 2- 



!ln(:/S) 



dlnR 



i?ISCi 

R 



(14) 



for this event, evaluated at i? = 10GM,/c^. I assume 
a pair of different initial surface density distributions for 
when the accretion flow is formed and shock-hcatcd. One 
is a (5-function annulus initially located at Ro = 0.05a w 
lO^GM./c^ from the SMBH. The second initial profile has 
a constant gas volume dcnsitj' (i.e., S oc -R"^), and extends 
from 0.5 -Ro to 1.52 Rq. The dimensions of the second pro- 
file was chosen so that both initial conditions share the same 
values for the total orbital angular momentum and gas mass. 

For both initial profiles, the accretion rate peaks at 
nearly the same time, and the subsequent decay behavior 
is nearly identical. For our viscosity prescription, the late- 
time mass accretion rate scales as M oc t~*^^. 



rflnM 
dint 



^d\nR 



(15) 



due to the late-time behavior of the central surface density 
profile as the fiow is depleted (Tanaka 2011). This is quite 
similar to the canonical late-time behavior M oc t^''^'^ for 
TDEs (Rees 1988; Phinney 1989), which arises from the or- 
bital energy distribution of the disrupted stellar material. As 
we will see below (§3.2), the light curve in any given emis- 
sion band does not follow any particular power-law (Lodato, 
King & Pringle 2009; GuiUochou & R amirez-Ruiz 2012; Shen 
& Matzner 2012, axgue the same for TDEs), and the emit- 
ted luminosity is of order the Eddington luminosity despite 
the accretion rate exceeding the critical value by an order of 
magnitude. 

A key model assumption with respect to the light curve 
is that the accretion flow forms promptly, i.e. on timescales 
shorter than ~ 10% of the binary orbital period. The fact 
that the dynamical times are much shorter where the stream 
shocks {R < 0.1a) favor this hypothesis, but this should be 
confirmed by numerical simulations. If the fueling accretion 
flow materializes on timescales longer than the fuel depletion 
timescale, then the light curves would rise and decay more 
gradually. 

In the bottom panel of Figure 2, 1 plot the fraction of the 
initial gas mass that is optically thick (r > 1). For flducial 
parameter values, the entire stream becomes optically thin 
on a timescale shorter than P, at which point the flow will 

have transitioned into a radiatively inefficient and under- 
luminous state. The amount of time it takes the gas to be 
depleted in this way depends on the total amount of gas 
in the flow. For very high accretion rates into the cavity, a 
quasi-steady state with small persistent discs around one or 
both SMBHs may result. 



3.2 Temperature and Energetics 

To model the emergent spectrum, I now estimate the radial 
temperature profile of the disc. I assume that the accreting 

gas is able to promptly dispose of the viscously generated 
heat via graybody radiation and advection: 




Figure 2. Top panel: The mass flux M evaluated aX R = 
10 GMtj(? as a function of time, for two initial density profiles: a 
(5 function centered at Ro = 0.05o (thick curves) and a constant- 
density annulus extending from 0.5i?o to 1.5i?o- The total mass 



of the stream is initially Mf^^ 



■ 0.05McrifP, where M^rit is the 



critical accretion rate for the secondary mass. Bottom panel: Tlie 
fraction of the initial mass in the accretion flow that is optically 
thick r > 1, assuming electron-scattering opacity. There are two 
(nearly overlapping) curves for each of the two initial conditions. 



where 



Qa,dv 



4 

: 2 X 



3r 



and 



M 
2^' 



12 



dluT 



dlnp 



dlnR dlnR 



(17) 
(18) 
(19) 



(16) 



represent viscous heating, radiative cooling and advective 
cooling (Abramowicz et al. 1988), respectively. Outer regions 
of the accretion flow (where the surface density is lower) that 
are optically thin are assumed to be radiatively inefficient 
and thus excluded from the emergent flux calculation. 

In this model, the viscosity is scaled to the gravitational 
potential via equations (7) and (9) as = 
While this approach has the advantage that the surface 
density can be evaluated easily by exactly solving a linear 
differential equation, it differs from the standard Shakura- 
Sunyacv prescription u = ac^/fiK in that the viscosity is not 
explicitly coupled to the thermal evolution of the gas (equa- 
tions 17 through 19). For the problem at hand, however, 
using the canonical viscosity prescription poses several prac- 
tical difficulties and is unlikely to lead to a better physical 
understanding of the flares. First, given the rapid evolution 
of the flares (as compared to the canonical viscous timescale, 
due to sharp gradients) , the wide range of densities and tem- 
peratures involved, and the inherent complexity of the flow 
geometry, the assumption of a constant effective a parame- 
ter is almost certain to be physically inaccurate. Second, the 
thin disc model is known to be canonically prone to thermal 
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and viscous instabilities (Lightman & Eardley 1974; Shakura 
& Sunyaev 1976; Pringle 1976; Piran 1978), which are active 
subjects of investigation in the context of the magnetohydro- 
dynamic studies of accretion flows (Hirose, Krolik & Blaes 
2009; Jiang, Davis & Stone 2013). If the accretion flow in 
question is thermally unstable, it should undergo runaway 
heating and be rapidly accreted (or blown away by its own 
radiative output), and thus the time evolution calculated 
here may be conservative. 

Given the above two theoretical difficulties and uncer- 
tainties associated with explicitly implementing the classical 
a prescription, the current approach is sufficient to demon- 
strate that the periodic streams can fuel rapid and luminous 
transient events. As a rough self-consistency check, I confirm 
that (i) the cooling/heating timescales of the accretion fiow 
are comparable to or longer than the timescale on which 
the entire accretion flow is depleted, and (ii) the evolved 
temperature profiles T{R, t), the disc aspect ratio H / R, and 
the derived kinematic viscosity a'csH are all consistent with 
equations (7) and (9) within a factor of a few throughout 
the accretion event. 

In figure 3, I present toy-model light curves, along with 
several snapshots of the emergent spectrum, for the accre- 
tion event modeled earlier in Figure 2. The binary mass is 
M,, = 10^ Mq, the period is P = 30 yr, and the accre- 
tion occurs onto the secondary SMBH with M, — 10^ Mq . I 
have included for reference the emission from the circumbi- 
nary disc, with m = 1, a = 0.1 and a cavity at i? < 2a. 

The top panel shows the bolometric light curve (black, 
solid curve) along with several monochromatic ones: R- 
band (700 nm; red, dotted curve), U-band (365 nm; blue, 
short-dashed curve), far UV (1500 A; magenta, long-dashed 
curve), and 0.5 keV (cyan, dash-dot curve). The thick lines 
show the light curves for the 5-function initial density pro- 
file, and the thing curves show those for the constant-density 
profile. Both light curves converge close to where the emis- 
sion peaks at roughly t = O.OIP. The emission drops off 
steeply as the gas is depleted and the outer accretion fiow be- 
comes optically thin and radiatively inefficient. None of the 
light curves follow a oc t"*'''^ power-law decay, even though 
the central accretion rate does (figure 2). The peak bolomet- 
ric luminosity is roughly equal to the Eddington luminosity 
of the secondary, even though the accretion rate is super- 
critical by an order of magnitude; this is consistent with 
steady-state models of slim accretion discs. 

In the bottom panel of figure 3, 1 show several snapshots 
of the SED, calculated as a graybody spectrum from the 
temperature profile as in the circumbinary SEDs in figure 1. 
On the lower-left corner is emission of the circumbinary disc, 
which in this example peaks in the infrared and is optically 
dim. This is the emission from the system most of the time, 
in between the flares. The composite (flare plus circumbi- 
nary disc) spectrum at the peak {t — O.OIP) and just as the 
flare is fading {t — 0.26P) are also shown to demonstrate 
that, in principle, the time evolution of the emission prop- 
erties can be quite dramatic. In addition to the viscous and 
thermal evolution of the accreting gas, the duration of the 
flare is also sensitive to the raw amount of mass deposited 
via the stream. 
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Figure 3. Toy light curves (top panels) and spectrum snapshots 
[bottom panels) for a M„ = 10^ Mq binary with a 1 : 10 mass 
ratio. The binary has a period of 30 yr, which places the peak 
emission frequency of the circumbinary disc in the infrared. In 
this example, the flares peak in the UV and soft X-ray frequencies 
and last several months before decaying rapidly over several years. 
Neither the bolometric nor the monochromatic light curves obey 
the t~^^^ decay of the central accretion rate. 



3.3 An alternate binary model for OJ287 

The BL Lacertae object OJ287 is known to produce optical 
outbursts every 11 — 12 years. The outbursts occur in pairs 
whose peaks are separated by approximately 2 years. 

Numerous studies have interpreted this system as a 
SMBH binary. In particular, Sillanpaa et al. (1988); Lehto & 
Valtonen (1996); Valtonen et al. (2006, 2008) have modeled 
the flares as resulting from a secondary SMBH impacting 
the accretion disc of the primary. In this disc-impact model, 
the secondary's orbit is elliptical and inclined with respect 
to the accretion disc, so that two impacts occur near peri- 
center. Valtonen, Ciprini & Lehto (2012) take the primary 
and secondary masses to be 1.8 x 10^° Mq and 1.2 x 10* Mq, 
respectively. 

Here, I interpret OJ 287 as an accreting binary with a 
period of 11.7 years (9 years in the binary's reference frame, 
for redshift z — 0.3) but with two substantive differences 
from the model of Valtonen et al. First, the binary's orbit 
is coplanar with a circumbinary accretion disc, and the disc 
has a circumbinary cavity with a radius of 2a. Second, I 
use a total binary mass of lO'' Mq, more consistent with 
independent observational results than the > lO^'' Mq value 
in the binary model (e.g., Xie et al. 2002; Liang & Liu 2003; 
Fan & Cao 2004; Gupta et al. 2012). The accretion rate in 
the circumbinary disc is taken to be rh = 5, by matching 
the optical-wavelength flux in between outbursts with the 
circumbinary disc output. (As noted above, because of the 
decreased flux from the center of the binary potential, the 
circumbinary disc can have m > 1 without being Eddington- 
limited.) I take the primary and secondary masses to be 8 x 
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Figure 4. A toy emission model for OJ 287. The system is mod- 



eled as a binary witli a total mass M, 



1O«M0, a 1 : 4 mass 



ratio, a redshift of z = 0.3 and an orbital period of P = 9 yr 
(J'obs = 11-7 yr). The binary's tidal torques are assumed to de- 
posit a stream near each SMBH, producing flares that are offset 
by 2yr (see text for details). Top panel: the bolometric (black) 
and U-band (blue) light curves, combing the emission from both 
flares and the circumbinary disc. Bottom panel: Snapshots of the 
SED during quiescence (marked with a vertical line and the letter 
"A" in the light curve), during the first flare ("B") and during 
the second flare ("C"). 



10 M0 and 2 x 10 M© , respectively. I suppose, arbitrarily, 
that they accrete 5% of m,, between them, and that the 
fraction consumed by each SMBH is inversely proportional 
to their mass ratio (the rest of the gas may contribute to an 
outflow, or cause a gradual pile-up and in the circumbinary 
gas). I take a' = 0.5 and Ro = 0.05a, and assume that the 
flares are offset by 2 years in the observer's frame, based on 
the observed delay between the two peaks. 

I present the light curves and SED snapshots of this 
model in flgure 4. The observable flux is computed from the 
source emission assuming a luminosity distance of 1600 Mpc 
to correspond to the source redshift. The bolometric lumi- 
nosity (black solid curves, top panel) fluctuates by an order 
of magnitude between quiescence and flaring. The U-band 
monochromatic flux (blue dotted curves) varies by two or- 
ders of magnitude. SEDs during quiescence ("A") and the 
peak of each flare ("B, C") are shown in the bottom panel. 
The flrst flare, due to accretion onto the secondary, lasts 
longer because the mass in the accretion flow is assumed to 
be larger. The steep dropoff in the SED of the circumbinary 
disc is similar to what is actually observed in OJ 287, but 
this feature could also be due to extinction by dust (Valto- 
nen, Ciprini & Lehto 2012). 

Many different spectra and light curves can be produced 
by varying the model parameters, but the production of two 
recurring flares at near-Eddington fluxes is a generic fea- 
ture for binaries with prograde accretion discs. Given the 
approximate nature of the model, I undertake no special ef- 



fort to obtain a detailed flt to the observed properties of OJ 
287. Rather, the goal of this exercise is merely to demon- 
strate that this simple picture — periodic streaming of gas in 
the disc cavity — can provide an explanation for the cadence, 
duration and energetics of the optical flares. With respect 
to OJ 287, the present model has several notable features. 
First, the interpretation does not require a primary mass of 
> 10^" ~ Mq , which is over an order of magnitude larger 
than several empirical estimates, as well as the M, — a rela- 
tion. Second, the cavity naturally explains the optical dim- 
ness of the system in between outbursts. In fact, the steep 
dropoff in the SED at optical frequencies arises simply from 
a rh ~ 1 disc around a M ~ 10^ M© binary being trun- 
cated at ii ~ 2a. Additionally, the initial formation of the 
flaring accretion flow by shock heating could explain the op- 
tical precursor flares observed before the major outbursts 
(Valtonen et al. 2006; Pihajoki et al. 2013). 

In the disc- impact model of Valtonen et al. (2006), the 
flares are caused by a precessing secondary impacting a cir- 
cumprimary disc that is stationary in space; the timing pre- 
dictions are strongly constrained by the binary parameters 
(the SMBH masses, spins and eccentricity give the preces- 
sion rate). In contrast, in the cavity-flare model I discuss 
here, some scatter in the flare timing is expected due to the 
intrinsic variability in the dynamics of the streams. Future 
observations could discriminate the two models by conflrm- 
ing or falsifying the binary orbit solution of the disc-impact 
model. 



4 POPULATION OF FLARING SMBH 
BINARIES 

How frequent are flaring events of this type? The forma- 
tion and dynamical evolution of sub-parsec SMBH binaries 
remain active topics of research in which guiding empirical 
facts are scarce (see, e.g., Colpi & Dotti 2011). However, 
one can make some educated estimates. For example, we 
know that a galaxy merger is a prerequisite for SMBH bi- 
nary formation, and we have a strong theoretical founda- 
tion for hierarchical structure formation. Fakhouri, Ma & 
Boylan-Kolchin (2010) found that the mean merger rates of 
dark matter halos in the Millennium Simulations (Springel 
et al. 2005; Boylan-Kolchin et al. 2009) can be described 
by a simple analytic fitting formula. For mergers with halo 
mass ratios between unity and 30 : 1, the merger rate (per 
halo per unit time) evaluates to 

0.13 



2.1 



1012 M, 



, dz 
It 



(20) 



per halo, for a given halo mass Mhaio and redshift z. I com- 
pute dz/dt in the standard ACDM cosmology with h = 0.70, 
= 0.72 and 0^ = 0.28 (Hinshaw et al. 2012). 
By relating the stellar velocity dispersion cr* of galaxies 
to the dark matter halo mass, one can attempt to convert 
the observed M, — a, relation (Ferrarese & Merritt 2000; 
Gebhardt et al. 2000; Giiltekin et al. 2009) to a relation be- 
tween M, and the host halo mass Mhaio (Ferrarese 2002; cf. 
Wyithe & Loeb 2003; Baes et al. 2003; Bandara, Crampton 
& Simard 2009), in a power-law form 
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(21) 
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Although this relation is known to have a large intrinsic scat- 
ter, it serves to provide an order-of-magnitude estimate. I 
choose {A, B, C) = (0.1, 5/3, 3/2), which are in rough agree- 
ment with the references above. Combining equations (20) 
and (21), one can estimate the rate at which two SMBHs of 
combined mass M,, end up in a single halo. 

Some fraction of such pairs will form binaries embed- 
ded in circumbinary discs. There are (semi-) analytic for- 
mulations of how such binaries evolve toward merger (Syer 
& Clarke 1995; Ivanov, Papaloizou & Polnarev 1999; Liu & 
Shapiro 2010; Haiman, Kocsis & Menou 2009). In the sim- 
plest description, the evolution can be written 



da 



: ttGW + Odis 



(22) 



At large separations, the binary evolves by depositing its 
orbital angular momentum into the circumbinary accretion 
disc (ddisc), and at close sepaxations by emitting gravita- 
tional waves (oiGw)- For the former, we are interested in 
secondary type-II migration whore the binary can open a 
central cavity. I take the simple prescription given by Syer 
& Clarke (1995) 



2rn(l + qf U 



(23) 



where g < 1 is the secondary-to-primary mass ratio, tv = 
2F? /{'iv) is the canonical viscous timescale evaluated at the 
cavity edge R — 2a and the dimensionless quantity k ~ 0.4 
depends on the surface density and viscosity profiles of the 
disc. I take v to be given by a standard Shakura-Sunyaev 
disc with a = 0.1 and m = 0.3 (typical for luminous AGN; 
e.g. KoUmeier et al. 2006), with the critical accretion rate 
evaluated with respect to the binary mass. I assume q — 0.1 
as a typical SMBH mass ratio in such binaries. The effects 
on the above prescriptions for varying the viscosity prescrip- 
tion and various system parameters (e.g., q, m, a) can be 
found in Haiman, Kocsis & Menou (2009); given the ap- 
proximate nature of this exercise, I will not explore this pa- 
rameter space. Note that Kocsis, Haiman & Loeb (2012a) 
find that the disc-driven migration of the binary can be sig- 
nificantly slower than the expression in equation (23); the 
actual lifetimes of SMBH binaries with circumbinary discs 
could be longer than estimated here. For the gravitational 
wave-driven evolution, I take 



ttGW 



64 



GM/c2 5 GM.. (l + qy \GM/c 



(24) 



from Peters (1964), assuming circular orbits for simplicity. 

The most liberal assumption that one can make is that 
once SMBH pairs are formed after a galaxy merger, all of 
them evolve toward coalescence while embedded in circumbi- 
nary discs. Thus, an upper limit for the average flaring rate 
per galaxy (Fflarc), averaged over all galaxies, can be esti- 
mated by multiplying the universal merger rate per halo 
with the total number of flares (i.e., orbits) produced by 
SMBH binaries as they evolve toward merger: 



-L merge {M..,z) 



da\ 
dt) 



da. (25) 



I've assumed one flare per binary orbit above, but it is trivial 
to generalize to two flares. Similarly, the majcimum fraction 



of galaxies that host a flaring SMBH binary between ob- 
served orbital periods Pi and P2 is 



/..(Pobs) 



e{M..,z) f 

J a 



a(f2,z) 



dt 



da. (26) 



The z-dependence of the above expressions assumes that the 
timescale for the formation of a disc-embedded binary takes 
place on timescales much shorter than a Hubble time. 

I conservatively choose the upper limit of integration 
in equations (25) and (26) to be aq = R{Q = l)/2, i.e. 
the binary semimajor axis where the Toomre parameter Q 
equals unity at the inner edge of the circumbinary disc. 
Outside R{Q = 1), the disc is thought to be gravitation- 
ally unstable and prone to fragmentation (e.g., Goodman 
2003). I also limit the integration to a minimum orbital pe- 
riod in the observer rest frame Pobs = (1 + z)P of 1 month. 
This cutoff is motivated by two considerations. First, Koc- 
sis, Haiman & Loeb (2012a) found that the cavity closes 
for rest-frame orbital periods shorter than ^ 10 days; this 
suggests that SMBH binaries would not produce flares on 
timescales shorter than weeks. Second, if a system were to 
produce flares on timescales shorter than 1 month, the peri- 
odic nature of the source should become transparent shortly 
after a single flare from the source is detected, i.e. dur- 
ing the course of standard monitoring of a TDE candidate. 
Here, I focus on flaxes whose recurrent nature would be chal- 
lenging to detect. Note that we can immediately place an 
order-of-magnitude estimate that /,, ~ 10~^, since we know 
Lmorgc ~ a few X Gyr^^ and that the timescale for a SMBH 
binary to evolve from oq to coalescence is ~ i,^ ~ Myr. 

I plot these upper limits for Fflare and /,, in figure 5, 
as a function of M,, and for several different subsets of bi- 
nary populations. From left to right, the vertical pairs of 
panels show binary populations at 2; = 0, z = 1 and z = 2. 
The top panels show (Fflare) (fiares per unit time, per all 
galaxies, in the galaxy rest-frame). The black, solid curves 
show all binaries with semimajor axes where Q = 1 at the 
cavity edge and where the observed orbital period equals 1 
month. The blue, dsished curves show the fiaring rate from 
only the binaries whose circumbinary discs are luminous in 
the V-bamd, and the red, dotted curves show those that are 
ultra-dim in the V-band in between flares. This was deter- 
mined by evaluating whether the redshifted peak frequency 
in the spectral emission of the circumbinary disc was higher 
or lower than 10^'' Hz (equation 5). Because a disc would 
emit emit much like a typical AGN below this cutoff and 
produce very little thermal emission above it, I do not ap- 
ply a luminosity cut. The figure shows that if all SMBH pairs 
merge via migration in a disc with a central cavity, then the 
event rate could exceed 10~^yr~^, compared to the TDE 
rate estimates of Ftde ~ IQ-'^yr-i (Ma gorrian & Tremaine 
1999). In other words, if ~ 1% of such pairs produced flares 
in the manner considered in this work, then the global event 
rate of such flares could be comparable to the rate of stel- 
lar TDEs. (However, recently merged systems may prefer- 
entially by obscured by dust.) The overall fraction of flares 
occurring in systems with optically dim circumbinary discs 
increases with z, as the peak disc frequency is redshifted to 
lower values. 

The bottom panels of figure 5 show the upper limits for 

the fraction of all galaxies that host an accreting SMBH 
binary in specified ranges of orbital periods Pobs in the 
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Figure 5. Top panels: Estimated upper limits of the flaring rate 
per galaxy, across all galaxies, in the source rest frame. The black, 
solid curves show the rate for all binaries with separations be- 
tween the value a gravitationally stable circumbinary disc can 
exist and the value where the observed binary period is 1 month. 
The blue, dashed curves show the rate only from binaries whose 
circumbinary discs are visible as optical AGN, and the red, dotted 
curves show the sources whose circumbinary discs are optically 
dim. These upper limits show that the global flaring rate could, in 
principle, exceed the rate of stellar TDEs. Bottom panels: Upper 
limits for the fraction of all galaxies that host flaring SMBHs. "All 
binaries" refers to the same range of semimajor axes as in the top 
panels, and the colored lines show the fractions for those binaries 
in the denoted ranges of orbital periods in the observer's rest- 
frame. For both the top and bottom panels, the redshift of the 
population increases from left to right as z = 0, 1, 2. 



observer's rest frame. Again, "all" binaries refers to those 
between where Q{2a) — 1 and Pobs = 1 rno. I have also 
plotted the maximum fraction of galaxies hosting SMBH 
binaries with periods 1 mo < Pobs < 1 yr (dotted red 
curves), 1 yr < Pobs < 10 yr (short-dashed blue curves), 
and 10 yr < Pobs < 100 yr (long-dashed purple curves). The 
figures suggest that flaring binaries should have orbital pe- 
riods (flare frequencies) of several months to several years 
if the binary mass is less than ~ lO^M©, and ^10 yr for 
M ~ W^Mq (coincident with the flare period of OJ 287). 

These upper limits are, of course, highly speculative, 
given the utter uncertainty of the formation and evolution of 
SMBH binaries, the fraction of these objects that form pro- 
grade circumbinary accretion discs, and the fraction of those 
objects that will be unobscured. These calculations do, how- 
ever, demonstrate the proof-of-concept that flaring SMBH 
binaries might be a relatively common occurrence, perhaps 
comparable to the incidence of luminous AGN activity. This 
provocative possibility should be falsifiable through obser- 
vations from upcoming surveys such as LSST or eROSITA. 



5 DISCUSSION 

Above, I have argued that an accreting SMBH binary 
with a prograde circumbinary disc may produce periodic 
UV/optical flares at cadences of years or decades, that this 
may occur in nuclei that are optically quiescent, and that 
such events may be common in the transient sky. That is 
good and well as a piece of educated speculation, but the 
inverse problem — that is, supposing that such flares are ob- 
served, then discerning whether the cause is an accreting 
SMBH binary — is a much more difficult one. 

It is well-known that most AGN exhibit aperiodic op- 
tical/UV fluctuations at the ~ 10% level over timescales of 
weeks to years (Ulrich, Maraschi & Urry 1997; Giveon et al. 
1999; Collier & Peterson 2001). However, examples in which 
the AGN appears to "turn on" briefly from a previously qui- 
escent state are rare (e.g., Gilli et al. 2000; Grupe, Thomas & 
Beuermann 2001). Transient AGN-like activity with rapidly 
decaying optical and UV light curves are generally thought 
to be stellar TDEs (Komossa & Greiner 1999; Komossa & 
Bade 1999; Gezari et al. 2006; Esquej et al. 2008; Gezari 
et al. 2009; Cappelluti et al. 2009; van Velzen et al. 2011; 
Cenko et al. 2012; Saxton et al. 2012). 

Indeed, the flares proposed in this work are physically 
similar to TDEs, in that an amount of gas comparable to 
that of a star is deposited close to a (possibly quiescent) 
SMBH — except that the "fuel" in TDEs is deposited at a 
distance of a few gravitational radii, whereas the gas streams 
leaked from the circumbinary disc are typically deposited at 
tens to thousands of gravitational radii (i.e., ^ 0.1a). One 
could also think of the flares as a gas-rich analogue of G2, 
the purported ~ 1O~^M0 gas cloud^ falling toward Sag A* 
(Gillessen et al. 2012). 

Disentangling TDEs from SMBH binaries is further 
complicated by the fact that the two phenomena are not 
mutually exclusive. On the contrary, it has been suggested 
that SMBH binaries — or a recently merged binary (Komossa 

6 Merritt 2008; Stone & Loeb 2011, 2012)— may even en- 
hance the rate of TDEs, without the aid of an accretion 
disc (Chen et al. 2009, 2011). Conversely, an accretion disc 
around a single SMBH may enhance the TDE rate with- 
out a binary (Karas & Subr 2007). While observing mul- 
tiple, TDE-like UV/optical transients in the same galactic 
nucleus in a human lifetime is thought to be extremely un- 
likely in the absence of a SMBH binary or merger remnant, 
an observation of multiple flares in a single nucleus (even if 
corroborated by additional circumstantial indicators such as 
morphological evidence of a recent galaxy merger) may not 
be sufficient to discern the nature of the source. 

Fortunately, the scenario proposed here has several ba- 
sic, predicted features (other than periodicity) that should 
help to distinguish flaring supermassive binaries from TDEs. 
The most signiflcant prediction is that the flares can be pro- 
duced by SMBHs with M. > IO^Mq, which are thought 
not to disrupt solar-type stars. Another is that it requires 
the presence of a circumbinary accretion disc with a cavity, 
which should have an anomalously red or optically dim SED 
that may be observed in the infrared; the cavity may also re- 
sult in unusually weak broad emission lines. AGN with such 



5 Murray-Clay & Loeb (2012) and Scoville & Burkert (2013) sug- 
gest that G2 is instead a low-mass star with an associated outflow. 
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properties are known to constitute a small fraction (~ 1%) 
of the general population (Gibson, Brandt & Schneider 2008; 
Shemmer et al. 2009). 

I have only discussed the thermal emission from the 
optically thick accretion flow. It is widely established that 
SMBH accretion can exhibit a whole zoo of activity. This di- 
versity is exhibited by TDE candidates, via the contribution 
of the outflowing wind (Strubbe & Quataert 2009, 2011), as 
well as the production of relativistic jets (Burrows et al. 
2011; Bloom et al. 2011; Zauderer et al. 2011) with asso- 
ciated radio transience (Giannios & Metzger 2011). Sesana 
et al. (2012) discussed the possibility that SMBH binaries 
may produce periodic variability in the Fe Ka emission line 
(but this may be diflicult to observe if the accretion flows 
inside the cavity are short-lived); the rapid onset of a strong 
accretion flow may trigger radiatively driven jets and out- 
bursts (Tanaka & Menou 2010). The accretion events dis- 
cussed in this paper, if they occur in nature, should have a 
wide variety of observable signatures. 

In constructing the toy model of the accretion flare, I 
made several ordcr-of-magnitude estimates, based on pre- 
vious studies, for quantities such as the mass deposited by 
each stream, and the distance from the SMBH where the 
stream circularizes and begins to accrete. Another key as- 
sumption is that the accretion flow is hot and optically 
thick. While the canonical thermal instability of radiation 
pressure-dominated a viscosity flows supports this assump- 
tion, the actual light curves and SEDs will depend on the 
detailed thermal and viscous evolution of the gas. The emis- 
sion signatures presented here arc based on idealized treat- 
ments of the accretion geometry and viscosity, and should 
be viewed only as a proof-of-concept demonstration. Finally, 
the toy model assumes that the fueling accretion flow forms 
on a short timescale much shorter than the binary orbital 
timescale; the peaks of the flares would be less punctuated 
if the fueling flows are slow to form and circularize. The va- 
lidity of the above theoretical treatments can be tested by 
future numerical simulations. 



6 CONCLUSION 

In this paper, I explored the possibility that the periodic 
leakage of circumbinary gas onto a SMBH binary exhibits 
luminous transient flares. The underlying mechanism behind 
this hypothesis consists of two characteristics of prograde 
accretion discs found in numerous simulations: (i) the pres- 
ence of a central, circumbinary cavity in a prograde accre- 
tion disc, which reduces the emergent flux of the system at 
frequencies above the optical; and (ii) the periodic, punc- 
tuated leakage of the disc gas into the cavity, which would 
trigger UV and soft X-ray emission as it accretes onto one 
or both SMBHs. The main flndings for these "flaring SMBH 
binaries" are as follows: 

(i) Given the orbital trajectories and accretion rates 
found in numerical simulations, whatever gas accretes onto 
the SMBH(s) will be depleted and become optically thin on 
a timescale shorter than the binary's orbital timescale, if the 
initial thermal energy of the accreting gas is comparable to 
its free-fall orbital energy. Thus, the streams may produce 
one or two flares per binary orbit, depending on whether 
one or both SMBHs develop optically thick accretion flows. 



The initial formation of the accretion flow via shock heating 

may also trigger a flare at lower frequencies. 

(ii) Unlike TDEs, these flares can be powered by SMBHs 
of arbitrary mass, including those in the lO^M© class. For 
low mass SMBH binaries with short orbital periods the ac- 
cretion cavity may close (Kocsis, Haiman & Loeb 2012a) 
and the recurring flares may not be produced. 

(iii) Because the presence of the cavity reduces the high- 
frequency thermal emission of the circumbinary disc, in be- 
tween the outbursts such a system may appear as an anoma- 
lously coloured AGN or be mistaken for a quiescent galac- 
tic nucleus. The circumbinary disc may also have unusually 
weak broad-line or soft X-ray emission. 

(iv) If the gas sound speed scales with its Keplerian ve- 
locity, then the corresponding a viscosity prescription will 

drive a central accretion rate that decays as M cx t~*^^. 

(v) Energetic considerations point to an optically thick, 
marginally geometrically thin, advective flow similar to a 
slim disc. A toy treatment suggests that the thermal emis- 
sion of the flare will generically peak in the UV or soft X- 
rays, with variability extending to optical frequencies. 

(vi) This can cause periodic, orders-of-magnitude optical 
fluctuations, such as those exhibited by OJ 287. This inter- 
pretation differs from the disc-impact model proposed by 
Valtonen et al. (2006) in several ways: the required SMBH 
mass is smaller ( 10^ Mq); the optical dimness in between 
outbursts is a direct consequence of the cavity; whereas the 
disc-impact model strongly constrains the flare timing based 
on the binary orbit solution, some scatter is expected in the 
cavity-flare model due to the accretion dynamics. Addition- 
ally, the initial shock that circularizes the accretion flow in 
the cavity-flare model could explain the "precursor flare" 
events observed in this system. 

(vii) Monochromatic and bolometric light curves are not 
expected to exhibit a simple power-law decay. 

(viii) Most flaring SMBH binaries will have periods of 
~ 0.1 — 100 yr, with more massive systems capable of more 

luminous outbursts having longer cadences. For binaries 
with M ^ lO^M©, the orbital periods can be months or 
shorter; however, as noted above such low-mass binaries may 
not produce flares because the cavity may close. 

(ix) Depending on the fraction of galaxy mergers that 
produce SMBH binaries with prograde circumbinary discs, 
the global rate of such flares may be quite high. Under the 
most liberal assumptions — that most galaxy mergers end 
with such binaries and that they are not preferentially ob- 
scured because of the recent merger — the global rate of flar- 
ing SMBH binaries may exceed the theoretical rate of TDEs. 



If they occur in nature, these flares should be ob- 
served serendipitously by future surveys such as LSST or 
eROSITA. In the meantime, additional constraints on the 
prevalence of this mechanism can be placed by monitoring 
(candidate) TDE host galaxies for recurrent activity. The 
results presented here reinforce the notion that accreting 
SMBH binaries may be a common source of astrophysical 
transients (Haiman, Kocsis & Menou 2009; Tanaka, Haiman 
& Menou 2010). 
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